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ABSTRACT 

We present the optical discovery and sub-arcsecond optical and X-ray localization of the afterglow of the 
short GRB 120804A, as well as optical, near-IR, and radio detections of its host galaxy. X-ray observations 
with Swift/XRT, Chandra, and XMM-Newton extending to 5t « 19 d reveal a single power law decline. The 
optical afterglow is faint and comparison to the X-ray flux indicates that GRB 120804A is "dark", with a rest- 
frame extinction of Ay ost ss 2.5 mag (at z ~ 1.3). The intrinsic neutral hydrogen column density inferred from 
the X-ray spectrum, Nn;mt(z = 1.3) w 2 x 10 22 cm" 2 , is commensurate with the large extinction. The host galaxy 
exhibits red optical/near-IR colors. Equally important, JVLA observations at « 0.9- lid reveal a constant flux 
density of g„(5.8GHz) = 35 ±4 fiJy and an optically-thin spectrum, unprecedented for GRB afterglows, but 
suggestive instead of emission from the host galaxy. The optical/near-IR and radio fluxes are well fit with the 
scaled spectral energy distribution of the local ultra-luminous infrared galaxy (ULIRG) Arp220 at z ~ 1.3, 
with a resulting star formation rate of ps 300 M Q yr" 1 . The inferred extinction and small projected offset 
(2.2 ± 1.2 kpc) are also consistent with the ULIRG scenario, as is the presence of a companion galaxy at a 
separation of about 1 1 kpc. The limits on radio afterglow emission, in conjunction with the observed X-ray and 
optical emission, require a circumburst density of n ~ 10~ 3 cm" 3 an isotropic-equivalent energy scale of g 7i i s0 m 
£k.iso ~ 7 x 10 51 erg, and a jet opening angle of 9j > 8°. The expected fraction of luminous infrared galaxies 
in the short GRB host sample is ~ 0.01-0.3 (for pure stellar mass and star formation weighting, respectively). 
Thus, the observed fraction of 2 events in about 25 hosts (GRBs 120804A and 100206A), provides additional 
support to our previous conclusion that short GRBs track both stellar mass and star formation activity. 
Subject headings: gamma rays: bursts 



1. INTRODUCTION 

Short-duration gamma-ray bursts (GRBs) occur in a wide 
range of environments that include elliptica l and star fo rm- 
ing galaxies in the field and in clusters (e.g., Berger 2011] and 
referenc es therein). These galaxies have redshifts of z ~ 0.1 
to > 1 dBergeretalJ 120071: lRowlinsonetal.1 120101) star for- 
mation rates of SF R < 0.01 to « 40 M yr" 1 (lBergeril2009l: 
IPerley et al.ll201 lal). and stella r masses ofM* ss 10 9 -4 x 10 11 
M ? dLeibler & Bergerll20i(jl) . These properties are sugges- 
tive of a progenitor population that tracks both stellar mass 
and star format ion activity (though w ith a significant delay 
of ~ 0.3 Gvn lLeibler & Bergerll2010l) . in agree ment with the 
popular compact object co alescence model (e.g., Eich ler et al.l 
Il98l lNaravari^tan[T992h . 

In a similar vein, short GRBs also exhibit a range of 
explosion properties, with isotro pic-equivalent energies of 
g-v.isn~gK.isn ~ 1Q 49 - 10 52 erg dBergerl 120071: \tM^2Mft 
Nysewand er et al.ll2009t) . j et opening angles of 9j ~ few to 
> 20 deg in a few cases dBurrows et al.l 120061 iGrupe et"ai1 
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[20061 ISoderberg et all 120061: IWatson et al.l l2006URm£etak 
| 2012|), and circumburs t dens ities of n < 1 cm" 3 dBerger et al. 
120051: ISoderberg et al.ll2006l iFong et al.ll2012l) . The distribu- 
tion of opening angles is of particular importance since it im- 
pacts the true energy scale and event rate. To date, all mea- 
surements or limits on 9j have relied on X-ray observations 
thanks to the relative brightness and high detection fraction of 
the afterglo w in the X-ray band compared to the o ptical and 
radio bands dNy sewander et al.ll2009t lBergeril20 1 Ol) . 

Here we present the optical discovery and sub-arcsecond 
localization of the optical and X-ray afterglow of the short 
GRB 120804A, as well as optical, near-IR, and radio detec- 
tions of its host galaxy. The afterglow data constrain the 
burst properties (gK,iso> Qj, «)■ The host galaxy observa- 
tions identify it as an ultra-luminous infrared galaxy (ULIRG) 
at a photometric redshift of z ~ 1.3, making GRB 120804A 
one of the most distant short bursts known to date. This 
is the first ULIRG host in the short GRB sample, exceed- 
ing even the lu minous infrared gala xy (LIRG) likely host of 
GRB 100206A dPerlev et al.l 1201 lal) . We present the after- 
glow and host galaxy observations in $2] extract the explo- 
sion properties in §|3] and determine the host galaxy photo- 
metric redshift and properties in $4] Throughout the paper 
we report magnitudes in the AB system (unless otherwise 
noted) , use a Galactic extinction value of E(B-V) 0.204 
mag fSchlafly & Finkbeinerl l201 ll) . and employ the standard 
cosmological parameters: Hq = 71 km s" 1 Mpc" 1 , Qa = 0.73, 
and rijvi = 0.27. 

2. OBSERVATIONS AND ANALYSIS 

GRB 120804A was discovered with the Swift Burst Alert 
Telescope (BAT) on 20 12 August 4 at 00:55:47.8 UT 
(Baumgart ner et all2012l) . and was also detected with Konus- 
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WIND dSakamotoetalJl2012l) . The burst duration is T 9l) = 
0.81 ±0.08 s (15-350 keV) with a fluence of F 7 = (8.8 ± 
0.5) x 10" 7 erg cm" 2 (15- 150 keV) and (1.45 ±0.30) x 10" 6 
erg cm" 2 (15- 1000 keV). A joint analysis of the BAT and 
Konus -WIND data indica tes a peak energy of E p = 135^9 
keV dSakamoto et al.ll2~012l) . The spectral lags are 16 ± 12 ms 
(15-25 to 50- 100 keV) an d -5 ±6 ms (25 - 5 to 100-350 
keV), typical of short GRBs dNorris et al.ll2012l) . 

Swift/X-ray Telescope (XRT) observations commenced 
about 78 s after the burst and led to the identification of a fad- 
ing source, located at RA=15 h 35 m 47.51 s , Dec=-28°46'56.9" 
with an unc ertainty of 1.4" radiu s (90% containment, UVOT- 
enhanced; lOsborne et al.l 12012b . Observations with the 
UV/Optical Telescope (UVOT) began about 97 s after the 
burst, but no counterpart was detected to a 3c r limit of > 21.4 
mag in the white filter (at St w 97 - 247 s; IChester & Lienl 
120121) . Optical and near-IR observations with GROND start- 
ing at St « 1 .5 hr also led to n on-detections with > 2 2 mag 
(griz) and > 20.6 mag (7-band: ISudilovskv et alj|2012l) . 

2.1. X-ray Observations 

We analyze the Swift/XRT data using the HEASOFT pack- 
age (v6.11) and latest calibration files with the standard fil- 
tering and screening criteria. We generate the 0.3- 10 keV 
count-rate light curv e following the procedure described in 
Marg utti et al.l (|2012), with a re-binning scheme that ensures 
a minimum signal-to-noise ratio of 4 for each temporal bin. 
The data comprise 33 s (St w 97- 130 s) in Windowed Tim- 
ing (WT) mode, and 18 ks in Photon Counting (PC) mode 
(St « 150-6 x 10 4 s). 

We fit the time-averaged WT spectrum with an absorbed 
power-law model (tbabs x ztbabs x pow in Xspec) using a 
Galactic neutral hydrogen column density of /Vh,mw ~ 9.3 x 
10 20 cm" 2 dKalberla et al.ll2005l) . The resulting spectral pho- 
ton index is T = (2.5 ± 0.3) and the excess neutral hydro- 
gen column density is Nn,mt = (3.5 ± 1.1) x 10 21 cm" 2 at 
Z = (C-stat = 92 for 122 degrees of freedom; uncertain- 
ties are la). From the time-averaged PC spectrum we infer 
T = (2.1 ±0.1) and jV H .int = (3.2±0.5)x 10 21 cm" 2 at z = 
(C-stat = 283 for 316 degrees of freedom). We adopt the 
latter value of /Ve.int in the time-resolved spectral analysis, de- 
signed to account for the source spectral evolution and the re - 
sulting count-to-flux conversion factor (Mar gutti et al.ll20l2h . 
The uncertainties arising from the flux calibration are properly 
propagated in the unabsorbed 0.3- 10 keV flux light curve. 
We note that at z ~ 1.3 (see ® the best-fit parameters are 
r w 1 .9 and the intrinsic neutral hydrogen column density is 
Nhm ~ 2 x 10 22 cm" 2 . 

We also analyze a Chandra ACIS-S observation obtained 
on 2012 August 13.45 UT (St « 9.41 d) wit h a total exposure 
time of 19.8 ks (PI: Troja: lTroja et al.ll2012l) . The X-ray after- 
glow is detected with a significance of about 13cr at a count 
rate of 1.6 ±0.3 cps in the 0.5-8 keV range. Adopting the 
best-fit spectral parameters from the XRT analysis this trans- 
lates to an unabsorbed 0.3- 10 keV flux of (2.8 ±0.5) x 10" 14 
erg s" 1 cm" 2 . 

We further obtained an XMM-Newton observation on 2012 
August 22.91 UT (St « 18.9 d; PI: Margutti) to search for the 
signature of a jet break at late time. We analyze the EPIC 
data with the XMM Science Analysis System (SAS vl 1.0.0), 
selecting events with PATTERN < 12 for the MOS cameras, 
PATTERN < 4 for the pn camera, and FLAG = for both. 
To reduce the contamination by soft proton flares, we screen 



the original event files using a sigma-clipping algorithm. The 
remaining good science time is 29.5 ks for MOS1 and MOS2, 
and 25.5 ks for pn. The X-ray afterglow is detected in the 
MOS1 and MOS2 images with (1.0±0.3) x 10" 3 cps (0.2- 
10 keV) and in the pn image at the with (4.3 ±0.7) x 10" 3 
cps (0.2- 10 keV). We perform spectral analysis using the 
evselect tool, with the response files generated with the 
rmf gen and arf gen tools. The resulting inter-calibration 
factor of MOS1 and MOS2 relative to pn is 0.9- 1.1. We 
adopt the value of A^H.int from the XRT analysis, leading to a 
pn unabsorbed 0.3- 10 keV flux of (2.4 ±0.5) x 10" I4 erg s" 1 
cm" 2 . 

2.2. Optical Afterglow Discovery and Relative X-ray 
Astrometry 

We obtained two epochs of /-band im aging with the Ge m- 
ini Multi-Object Spectrograph (GMOS; Hook et al.ll2004l) on 
the Gemini-North 8-m telescope on 2012 August 4.27 and 
7.30 UT (St w 0.23 d and » 3.26 d, respectively). The ob- 
servations consisted of 1980 s and 2880 s, respectively, in 
0.65" seeing. We process the data using the gemini pack- 
age in IRAF, and perform photometry using a zero-point of 
28.46 ±0.10 mag measured on the nights of August 2-8 UT. 
We perform digital i mage subt raction of the two epochs with 
the ISIS package (Alard 2000), and recover a fading source 
with m, = 26.2 ±0.2 AB mag (with an additional system- 
atic uncertainty of ±0.1 mag due to the zero-point uncer- 
tainty). We consider this source to be the optical afterglow 
of GRB 120804A. Corrected for Galactic extinction the re- 
sulting flux density is F u = 0. 17 ± 0.05 /iJy. Images of the two 
epochs and the resulting subtraction are shown in Figure Q] 

We determine the absolute position of the afterglow by 
astro metrically matching the images to the 2MASS refer- 
ence frame using 45 common sources. The resulting astro- 
metric uncertainty is 0.15" in each coordinate. The after- 
glow position in the residual image is RA=15 h 35 m 47.479 s , 
Dec=-28°46'56.17" (J2000), with a centroid uncertainty of 
about 0.05" in each coordinate. 

To locate the Chandra X-ray afterglow on the optical im- 
ages we perform differential astrometry using 4 common 
sources. We find a relative offset between the two coordi- 
nate frames (Chandra to Gemini) of <5RA = -0.06 ±0.21" and 
(5Dec = -0.01 ±0.17", leading to a refined X-ray afterglow 
position of RA=15 h 35 m 47.478 s , Dec=-28°46'56.30" with an 
uncertainty of about 0.35" radius that takes into account a cen- 
troid uncertainty of about 0.06"; see Figure Q]). The Chandra 
position is in excellent agreement with the optical afterglow 
position. 

2.3. Optical/Near-IR Host Galaxy Observations 

In the second Gemini observation we detect an ex- 
tended source near the afterglow position with m, = 24.80 ± 
0.15 mag (including zero-point uncertainty and corrected 
for Galactic extinction) located at RA=15 h 35 m 47.477 s , 
Dec=-28°46'56.44", with a centroid uncertainty of about 
0.10". We consider this source to be the host galaxy of 
GRB 120804A. Given the brightness of the galaxy, the prob- 
ability of chance coincidence using a radius of 1" (e.g., 
iBloom et alJl2002l;lBergerl2010l) is P cc 0.02. We investigate 
potential association with brighter galaxies (« 19-21 mag) at 
larger offsets (w 0.2- 1') but find chance coincidence proba- 
bilities of P cc w 0.2-0.7, indicating that these are not likely 
to be associated with GRB 120804A. 
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We also observed GRB 120804A with GMOS on the 
Gemini-North telescope on 2012 August 4.24 UT (St « 0.20 
d) in r-band with a total exposure time of 2340 s. We pro- 
cess the data using the gemini package in IRAF, and per- 
form photometry using a zero-point of 28.41 ±0.02 mag mea- 
sured on the nights of July 27 - August 9 UT. Photome- 
try at the position of the host galaxy reveals a faint source 
with m r = 26.2 ± 0.2 mag, corresponding to a flux density of 
F v = 0.19 ±0.04 /iJy (corrected for Galactic extinction). Since 
this was an early observation, the measurement may be con- 
taminated by afterglow emission. However, taking into ac- 
count a spectral shape of F v oc v~ i) (l and the large rest-frame 
extinction (see © we conclude that the afterglow contribu- 
tion is sub-dominant, « 0.03 //Jy. 

We observed the host galaxy in /-band with the FourStar 
near-IR camera on the Magellan/Baade 6.5-m telescope on 
2012 August 28.98 UT with a total on-source time of 2390 s. 
We analyze the data using a custom pipeline in python, and 
perform photometry using common sources with the 2MASS 
catalog. The host galaxy has a measured brightness of my = 
23.05 ± 0.20 mag (corrected for Galactic extinction). 

We obtained Y- and /T,-band observations with the High 
Acuity Wide-field K-band Imager (HAWK-I) on the VLT 
starting on 2012 September 7.96 UT, with a total exposure 
time of 1320 s in each filter. We produce dark-subtracted 
and flat-fielded images of the field using the HAWK-I pipeline 
within esorex, and perform photometry on the /^-band im- 
age relative to the 2MASS catalog. We determine the T-band 
calibration using the instrumental zero-point (appropriate for 
our observations, which were obtained in good conditions), 
and confirm that this is appropriate by extrapolating 2MASS 
/-band photometry to observations taken in the /-band. The 
host galaxy is detected in the /T r band with thk = 22.0 ±0.1 
mag, and weakly in T-band with my = 23.7 ± 0.3 mag (both 
values are corrected for Galactic extinction). 

Finally, we obtained spectroscopic observations of the host 
galaxy with the FOcal Reducer and low dispersion Spectro- 
graph (FORS2) on the VLT on 2012 August 19.01 UT. The 
observations consisted of 4 x 600 s exposures, covering the 
wavelength range 4300 - 9300 A. We set the slit position an- 
gle to 149° to cover the host and nearby galaxy (see Figure[T|i. 
We detect no continuum or line emission at the position of 
the host galaxy, and only a faint continuum from the nearby 
galaxy. 

The host galaxy photometry is summarized in Table [2] We 
also provide photometry for the extended source located w 
1.4" to the south-east of the host galaxy position (Figure[TJ. 

2.4. Radio Observations 

We observed GRB 120804A with the Karl G. Jansky Very 
Large Array (JVLA) starting on 2012 August 4.97 UT (St « 
0.93 d) at a mean frequency of 5.8 G Hz. We utilized the 
WIDAR correlator dPerlev et alj|201 lbl) with a bandwidth of 
about 1 GHz in each sideband, centered at 4.9 and 6.7 GHz. 
All observations were undertaken in the B configuration, uti- 
lizing 3C286 for bandpass and flux calibration, with inter- 
leaved observations of J 1522-2730 for gain calibration. We 
calibrate and analyze the data using standard proced ures in 
the A stronomical Image Processing System (ALPS; iGreisenl 
2003), and list the resulting flux density measurements in Ta- 
bled] 

In the first four epochs (St rj 0.9- 10.9 d) we detect a sin- 
gle unresolved source coincident with the optical and X-ray 
afterglow, as well as the host galaxy: RA=15 h 35 m 47.485 s 



(±0.008), Dec=-28°46'56.44" (±0.20). However, the 
source remains constant in brightness, with flux densities of 
F u (4.9 GHz) = 43 ± 4 ^Jy and ^(6.7 GHz) = 25 ± 4 /jjy. The 
steady brightness and optically thin spectrum, F v oc v~ x ,7±0 ,8 , 
are unprecedented fo r GRB radio afterglows at early time 
dGranot & S ari 2002), but are instead suggestive of emission 
from the host galaxy. Indeed, the centroid of the radio source 
shows a smaller offset relative to the optical host galaxy cen- 
troid (0.10") than to the optical afterglow position (0.28"). 
The probability of chance coincidence for a source of this 
brightness within ~ 1" of the optical/near-IR host galaxy po- 
sition is P cc m 2 x IP" 4 , based on the number counts of faint 
5 GHz sources (Fomalont et al.lfl991h . This indicates that the 
radio source is the host galaxy of GRB 120804A. The obser- 
vation on 2012 September 1 1 UT was obtained during a JVLA 
re-configuration leading to poorer noise characteristics in the 
resulting image. The 3a limits are consistent with the earlier 
detections of the steady source. 

We use the lack of variability to place a 3a upper limit on 
the radio afterglow brightness of /v(5.8GHz) < 20 /iJy. 

In addition, we obtained JV LA observations of the LIRG 
host galaxy of GRB 100206 A (Pe rlev et al.ll20TTah on 2012 
September 10.27 UT to determine whether its large total star 
formation rate inferred from optical/IR data (~ 40 M© yr" 1 ) 
produces radio emission. The observing setup and data analy- 
sis follow the procedure described above. We observed 3C48 
for flux and bandpass calibration, and interleaved observa- 
tions of J0238+1636 for gain calibration. There is a bright, 
contaminating source in the field (27 mJy at 4.9 GHz and 
7 mJy at 6.7 GHz), located about 4.5' from the position of 
GRB 100206 A. We therefore image the field utilizing self- 
calibration techniques on this bright source, leading to a non- 
detection of radio emission with a conservative 5a limit of 
F„(5.8GHz) <0.1 mJy 

3. AFTERGLOW PROPERTIES 

We model the X-ray data, optical detection, and radio up- 
per limits using the standard afterglow synchrotron model 
dGranot & Saril20 02). We follow the standard assumptions of 
synchrotron emission from a power-law distribution of elec- 
trons (iV(7) oc j~ p for 7 > 7„,) with constant fractions of the 
post-shock energy density imparted to the electrons (e e ) and 
magnetic fields (e#). The additional free parameters of the 
model are the isotropic-equivalent blast-wave kinetic energy 
CEic,iso)> an d the circumburst density (n for a constant density 
medium). 

We note three key observational facts to guide the after- 
glow modeling. First, the X-ray light curve is best fit with a 
single decline rate of ax = -0.93 ± 0.06, which coupled with 
the spectral index of (3 X ~ -0.9 ± 0. 1 (ED indicates that the 
synchrotron cooling frequency is v c ~ vx and that p w 2.2. 
Second, the peak of the X-ray afterglow is ss 10-20 p,iy at 
St « 200 s. This suggests that the radio light curve will even- 
tually reach a similar peak flux density if the synchrotron peak 
frequency is v m « v x at St sa 200 s. Such a peak flux den- 
sity is consistent with our inferred radio upper limits. Finally, 
the optical and X-ray flux densities at St « 5.5 hr are com- 
parable, indicating that the optical to X-ray spectral index is 
Pox ~ 0, compared to an expected slope of f» 0.6 (for p = 2.2 
and v c ~ vx). This shallow slope indicates that GRB 120804A 
can be classified as a " dark" burst, following the definition of 
Uakobsson etail (120041) . 

Guided by these results, we find that the data can be fit 
with the following parameters (using z - 1.3; see ©: £k,iso ~ 
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8 x 10 51 erg, n ps 10" 3 cm -3 , e c ps 0.3, e B ps 0.1, and p ps 2.2. 
The inferred blast-wave energy is comparable to the isotropic- 
equivalent 7-ray energy, E 7 j so ps 6 x 10 51 erg. To explain the 
suppressed optical emission we also require A{j? st ps 2.5 mag 
(at z = 1.3). The resulting light curves are shown in Figure[2] 
Models with higher values of Ekmo and/or n lead to a lower 
value of v m and a larger peak flux density, and therefore vio- 
late the radio limits (see for example dotted line in Figure|2). 

In addition, the lack of a break in the X-ray light curve to 
at least ps 20 d places a lower bound on the jet collimation 
angle, with Oj > 8°, where we use the values of £k,iso and 
n inferred above. This indicates a beaming correction factor 
of 1 = [l-cos(6'j)]~ I < 120, and hence a beaming-corrected 
energy scale of E y +E K ps 1 . 1 x 10 50 - 1 .4 x 10 52 erg; the upper 
bound is set by isotropy. 

Since in the afterglow model above, the X-ray light curve 
sets the overall flux density scale, an alternative explanation 
for the radio non-detections and the low optical flux den- 
sity is that the X-ray emission is dominated by a different 
emission mechanism. One possibility is co ntribution from 
inverse Compton emission (Sari&Esin 200l|), but this re- 
quires a large density of n > 10 2 cm" 3 . Another possibility is 
emission from a newly-b orn rapidly-spinning magnetar (e.g., 
Zhang & Meszaros 2001), but the expected evolution in this 
case is a relatively constant brightness for a duration similar 
to the spin-down timescale, followed by F v oc f" 2 at later time 
(for a typical braking index of 3). This is quite distinct from 
the observed single power law of F v oc f" 1 at St > 200 s. 



4. A ULIRG HOST GALAXY 

The observed optical/near-IR spectral energy distribution 
(SED) of the host galaxy exhibits a red color of i — K ps 2.8 
mag (ps 4.3 Vega mag); the less certain r — K color is ps 3.9 
mag (ps 5.6 Vega mag). There is also noticable steepening 
between the i- and F-band filters. These properties are in- 
dicative of a Balmer/4000 A break at z ~ 1.3 (with a range of 
ps 1.1-1.5) and either an old or dusty stellar population. 

Taken in conjunction with the radio detection, the SED is 
reminiscent of ULIRGs (Lhr > 10 12 L Q ). To investigate this 
possibility we compare the host galaxy fluxes to the SED 9 of 
the local ULIRG Arp220. As shown in Figure [3j at z *» 1.3 
a simple scaling of Arp 220 provides a remarkable fit to the 
data. An elliptical galaxy template (2 Gyr old population; 
Pollet taet alJ l2007h at z ps 1.3 provides a reasonable fit in 
the optical/near-IR (although it underestimates the observed 
r- and A^-band fluxes), but cannot explain the radio emission 
(Figure [3]). 

At z ~ 1.3 the rest-frame B-band absolute magnitude is 
Mb ~ -20.2 mag, or Lg ps 0.2 L* in comparison to the B-band 
luminosity function at z ~ 1 .3-2 dllbert et al.l2005l) . The rest- 
frame /f-band absolute m agnitude is Mk ps -2 2. 1 mag (-24.0 
Vega mag, or L K ps 0.4L*: ICaputi et al.l2006bl) . This luminos- 
ity corresponds to a stellar mass of M*5 x ps 10 10 (using a 
chara cteristic mass-to-light ratio of M*/Lk ps 0.3; lDrorv et al.l 
2004). The infrared bolometric luminosity scaled using the 
SED of Arp 220 is L FIR ps 10 12 L . 

The unobscured star formation rate is inferred from the ob- 
served r-band (Aq ~ 2700 A) to be ps 1 M Q yr" 1 . However, 
from the observed radio emission we estimate the total star 



formation rate to be much larger dYun & Carillill200"2h : 

F u,^ydl tQpc /(l+Z) 



SFR : 



2 < ,,-0.75 , n 7 ,. -0.1 



■ 300 Mq yr" 



(1) 



9 Obtained from the SWIRE Template Library iPolletta et al. 2007j). We 
note that ULIRG templates exhibit some variability, but this will have little 
effect on the resulting photometric redshift. 



where vq = (1 +z)f bs is the rest frequency, and di is the lumi- 
nosity distance. Assuming that the host stellar mass has been 
assembled with this star formation rate gives a characteristic 
age of about 0.15 Gyr. Using the limit on radio emission from 
GRB 100206A in Equation[f]we find SFR < 60 M Q yr" 1 . 

The large star formation rate in the host of GRB 120804A 
is also expected to pr oduce X-ray emission, with Lx ps 7 x 
10 39 x SFR erg s" 1 dWatson et all 12004 IVattakunnel et al l 
|2012|) . For the values inferred above we find an expected lu- 
minosity of Lx ps 2 x 10 42 erg s" 1 , corresponding to a flux 
of Fx ps 4.5 x 10~ 16 erg s" 1 cm" 2 . This is about 55 times 
lower than the afterglow flux measured with XMM-Newton 
at St ps 19 d, and hence consistent with a star formation origin 
for the radio emission. 

An alternative interpretation of the radio detection is emis- 
sion from an active galactic nucleus (A GN). Matching the 
standard template of a radio-quiet AGN (Shang et alj|201 lh 
to the observed radio flux density (using z = 1.3), leads to 
a substantial over-estimate of the optical/near-IR brightness, 
and a much bluer color (Figure |3). Similarly, the expected 
X-ray flux in this scenario is Fx ps 6.4 x 10" 14 erg cm" 2 s" 1 , 
several times brighter than the XMM-Newton detection of 
the afterglow. Thus, we can rule out a radio-quiet AGN ori- 
in. If we instead use the standard radio-loud AGN template 
Shang ^tai1l201 ll) matched to the radio flux density, we find 
that the expected optical/near-IR and X-ray fluxes are ~ 1 — 2 
orders of magnitude fainter than measured (Figure[3]l. Indeed, 
the observed lower bound on the ratio of radio to X-ray lumi- 
nosity, vL vm i/Lx >2x 10~ 4 , along with an upper bound of 
Li < 5 x 10 44 erg s" 1 are in good agreement with sam ples 
of low-luminosity AGN (e.g.,[Terashima & Wilson 20Q3j). In 
this scenario, the optical/near-IR emission is instead domi- 
nated by a stellar component, either from an old population 
or a reddened young population with a modest star formation 
rate of a few M Q yr" 1 . 

While both a ULIRG and an AGN origin can explain the 
observed radio emission, we note that the former also offers 
a natural explanation for the inferred extinction, and fits the 
broad-band host SED with a single component. We therefore 
consider a ULIRG as the more likely explanation for the host 
galaxy of GRB 120804A. 

The offset between the optical afterglow position and host 
galaxy centroid is 0.27 ±0.15", corresponding to 2.2 ± 1.2 
kpc at z ps 1.3. This is a relativ ely small offset for short GRB s 
(with a median of ps 5 kpc; iFong et alJl20Tot lBergerll20ldl) , 
though not unprecedented. Still, in the context of a ULIRG 
origin, the small offset is consistent with the inferred after- 
glow extinction. 

Finally, we determine a photometric redshift for the galaxy 
located ps 1 .4" from the host galaxy position and find z ps 1.3, 
consistent with that of the host galaxy (Figure [3). Thus, it 
appears likely that these galaxies, with a projected separation 
of about 1 1 kpc, are interacting or merging. This is not un- 
expected since at least some ULIRG activity is triggered by 
galaxy mergers. 

5. SUMMARY AND CONCLUSIONS 

We presented the discovery and sub-arcsecond localization 
of the optical and X-ray afterglow of the short GRB 120804A. 
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A comparison of the observed fluxes points to substantial rest- 
frame extinction of Ay° st rs 2.5 mag, commensurate with the 
large neutral hydrogen column density, Na,mt ~ 2 x 1 22 cm -3 . 
In conjunction with deep radio limits we infer an energy of 
fiUiso ~ £k.iso ~ 7 x 10 51 erg and a low circumburst density 
of n ~ 10" 3 cm" 3 . The lack of a break in the X-ray afterglow 
at < 20 d, leads to 9j > 8°, in line wit h existing measure- 
ments of short GRB jets ( iFong et al.ll2012l) . We note that a 60 
ks Chandra observation obtained at St rj 46 d (PI: Burrows) 
may extend this limit to > 1 1° if no break is detected (i.e., an 
expected flux of rs 7 x 10~ 15 erg s" 1 cm" 2 ), or determine an 
angle of 9j ri 10° if a break at an intermediate time of rs 33 
d is detected (i.e., leading to an expected flux of rj 5 x 10~ 15 
erg s" 1 cm" 2 ). In either case, our estimate of 9j > 8° will not 
change significantly. 

We also detect the host galaxy of GRB 120804A, which ex- 
hibits red optical/near-IR colors and radio emission that are 
well-matched by the SED of a ULIRG (Arp220) at z w 1.3. 
The inferred total star formation rate is ks 300 M Q yr" 1 . A 
low-luminosity radio-loud AGN in an elliptical galaxy can- 
not be definitively ruled out, but the ULIRG interpretation, 
combined with the small projected offset of 2.2 ±1.2 kpc, 
more naturally explains the inferred afterglow extinction. The 
host galaxy is part of an interacting/merging system, which 
is not unexpected for ULIRGs. Observations with the Ata- 
cama Large Millimeter/submillimeter Array (ALMA) will be 
able to robustly distinguish the two scenarios. At 240 GHz 
(ALMA band 6), the expected flux densities are about 0.3 mJy 
and 1 /iJy, for the ULIRG 10 and AGN cases, respectively. The 
former can be detected with high significance in a short obser- 
vation. Similarly, the 60 ks Chandra observation at St « 46 d 
can achieve a limiting flux of Fx ~ 1.5 x 10~ 15 erg s" 1 cm -2 , 
which is still about 3 times higher than the expected X-ray 
emission due to star formation. Therefore, if flattening of the 
X-ray light curve is detected in this observation, it may point 
to an AGN contribution. 

The host galaxy of the short GRB 100206A at z = 0.407 
is also a luminous infrared galaxy, L IR r* 4 x 10 11 L Q , with 
SFR w 30-40 Mfp yr" 1 inferred fro m SED fitting of data at 
- 0.3-10 fim dPerlev etalj|2011al) . Here we find SFR < 60 
M© yr" 1 from radio observations. Of the ri 25 short GRBs 
with robust host galaxy associations (Fong et al. in prepara- 
tion), the hosts of GRBs 120804A and 100206 A are the only 
galaxies with clear LIRG/ULIRG properties. At z ~ 1, the 
population of ULIRGs and bright LIRGs accounts for ~ 25% 

10 This flux density is relevant for the scaled SED of Arp 220. Variations 
in the dust temperature would lead to a range of potential flux densities (e.g., 



of the total star for mation rate density dLe Floc'h et al.ll2~005l: 
Caputi et al. 2007), but their contr ibution to the stella r mass 
density is small, ~ 1 -few% (e.g. JCaputi et al. 2006a). Thus, 
the fraction of ~ 5- 10% in the short GRB host population 
is intermediate between the star formation and mass weighted 
fractions. This is expected for a progenitor population that 
tracks both star formation and stel lar mass, consistent w ith 
previous findings for short GRBs (Leibler & Berger 2010). 



We thank Karina Caputi, Ranga Chary, Francesca Civano, 
Martin Elvis, and Jane Rigby for useful discussions about 
ULIRGs and AGN. We also thank Brian Metzger for infor- 
mation on the magnetar scenario. The Berger GRB group 
at Harvard is supported by the National Science Foundation 
under Grant AST-1 107973, and by NASA/Swift A07 grant 
NNX12AD69G. E.B. acknowledges partial support of this re- 
search while in residence at the Kavli Institute for Theoretical 
Physics under National Science Foundation Grant PHY11- 
25915. The Dark Cosmology Centre is funded by the Danish 
National Research Foundation. V. M. acknowledges funding 
through contract ASI-INAF 1/004/11/0. Observations were 
obtained with the JVLA (program 12A-394) operated by the 
National Radio Astronomy Observatory, a facility of the Na- 
tional Science Foundation operated under cooperative agree- 
ment by Associated Universities, Inc. The paper includes 
data gathered with the 6.5 meter Magellan Telescopes lo- 
cated at Las Campanas Observatory, Chile; with ESO Tele- 
scopes at the La Silla Paranal Observatory under programme 
ID 089.D-0450; with the Gemini Observatory, which is op- 
erated by the Association of Universities for Research in As- 
tronomy, Inc., under a cooperative agreement with the NSF on 
behalf of the Gemini partnership: the National Science Foun- 
dation (United States), the Science and Technology Facilities 
Council (United Kingdom), the National Research Council 
(Canada), CONICYT (Chile), the Australian Research Coun- 
cil (Australia), Ministerio da Ciencia, Tecnologia e Inovacao 
(Brazil) and Ministerio de Ciencia, Tecnologia e Innovation 
Productiva (Argentina). This work is based in part on observa- 
tions obtained with XMM-Newton, an ESA science mission 
with instruments and contributions directly funded by ESA 
Member States and the USA (NASA). 

Facilities: Swift (XRT), CXO (ACIS-S), XMM-Newton 
GeminkNorth (GMOS), Magellan (FourStar), VLT (HAWK- 
I), JVLA 

IYun&Cariirl2 002; Mic halowski et aU2008l) . 



REFERENCES 



Alard, C. 2000, A&AS, 144, 363 

Baumgartner, W. H., et al. 2012, GRB Coordinates Network, 13581, 1 

Berger, E. 2007, ApJ, 670, 1254 

Berger, E. 2009, ApJ, 690, 231 

Berger, E. 2010, ApJ, 722, 1946 

Berger, E. 2011, New Astronomy Reviews, 55, 1 

Berger, E., et al. 2007, ApJ, 664, 1000 

Berger, E., et al. 2005, Nature, 438, 988 

Bloom, J. S., Kulkarni, S. R., & Djorgovski, S. G. 2002, AJ, 123, 1111 
Burrows, D. N, et al. 2006, ApJ, 653, 468 

Caputi, K. I., Dole, H., Lagache, G., McLure, R. J., Dunlop, J. S., Puget, J.-L., 

Le Floc'h, E., & Perez-Gonzalez, P. G. 2006a, A&A, 454, 143 
Caputi, K. I., et al. 2007, ApJ, 660, 97 

Caputi, K. I., McLure, R. J., Dunlop, J. S., Cirasuolo, M., & Schael, A. M. 

2006b, MNRAS, 366, 609 
Chester, M. M., & Lien, A. Y. 2012, GRB Coordinates Network, 13582, 1 



Drory, N., Bender, R., Feulner, G., Hopp, U., Maraston, C, Snigula, J., & 

Hill, G. J. 2004, ApJ, 608, 742 
Eichler, D., Livio, M., Piran, X, & Schramm, D. N. 1989, Nature, 340, 126 
Fomalont, E. B., Windhorst, R. A., Kristian, J. A., & Kellerman, K. I. 1991, 

AJ, 102, 1258 
Fong, W., Berger, E., & Fox, D. B. 2010, ApJ, 708, 9 
Fong, W., et al. 2012, ApJ, 756, 189 
Granot, J., & Sari, R. 2002, ApJ, 568, 820 

Greisen, E. W. 2003, Information Handling in Astronomy - Historical Vistas, 
285, 109 

Grape, D., Burrows, D. N., Patel, S. K., Kouveliotou, C, Zhang, B., 
Meszaros, P., Wijers, R. A. M., & Gehrels, N. 2006, ApJ, 653, 462 

Hook, I. M., J0rgensen, I., Allington-Smith, J. R., Davies, R. L., Metcalfe, 
N., Murowinski, R. G., & Crampton, D. 2004, PASP, 116, 425 

Ilbert, O., et al. 2005, A&A, 439, 863 



6 



Jakobsson, P., Hjorth, J., Fynbo, J. P. U., Watson, D., Pedersen, K., Bjornsson, 

G., & Gorosabel, J. 2004, ApJ, 617, L21 
Kalberla, P. M. W., Burton, W. B., Hartmann, D., Arnal, E. M., Bajaja, E., 

Moras, R., & Poppel, W. G. L. 2005, A&A, 440, 775 
Le Floc'h, E., et al. 2005, ApJ, 632, 169 
Leibler, C. N., & Berger, E. 2010, ApJ, 725, 1202 

Margutti, R., Zaninoni, E., Bernardini, M. G., Chincarini, G., & for the Swift- 

XRT team. 2012, ArXiv e-prints 
Michalowski, M. J., Hjorth, J., Castro Ceron, J. M., & Watson, D. 2008, ApJ, 

672, 817 

Nakar, E. 2007, Phys. Rep., 442, 166 

Narayan, R., Paczynski, B., & Piran, T. 1992, ApJ, 395, L83 
Norris, J., Gehrels, N., & Barthelmy, S. D. 2012, GRB Coordinates Network, 
13585, 1 

Nysewander, M, Fruchter, A. S., & Pe'er, A. 2009, ApJ, 701, 824 
Osborne, J. P., Beardmore, A. P., Evans, P. A., & Goad, M. R. 2012, GRB 

Coordinates Network, 13577, 1 
Perley, D. A., Modjaz, M., Morgan, A. N., Cenko, S. B., Bloom, J. S., Butler, 

N. R., Filippenko, A. V., & Miller, A. A. 201 la, ArXiv e-prints 
Perley, R. A., Chandler, C. J., Butler, B. J., & Wrobel, J. M. 201 lb, ApJ, 739, 

LI 



Polletta, M., et al. 2007, ApJ, 663, 81 

Rowlinson, A., et al. 2010, MNRAS, 408, 383 

Sakamoto, T., et al. 2012, GRB Coordinates Network, 13614, 1 

Sari, R., & Esin, A. A. 2001, ApJ, 548, 787 

Schlafly, E. F, & Finkbeiner, D. P. 2011, ApJ, 737, 103 

Shang, Z., et al. 2011, ApJS, 196, 2 

Soderberg, A. M., et al. 2006, ApJ, 650, 261 

Sudilovsky, V., Elliott, J., & Greiner, J. 2012, GRB Coordinates Network, 
13579, 1 

Terashima, Y., & Wilson, A. S. 2003, ApJ, 583, 145 

Troja, E., Sakamoto, T., Gehrels, N, Donato, D., & Racusin, J. L. 2012, GRB 

Coordinates Network, 13640, 1 
Vattakunnel, S., et al. 2012, MNRAS, 420, 2190 

Watson, D., Hjorth, J., Jakobsson, P., Pedersen, K., Patel, S., & Kouveliotou, 

C. 2004, A&A, 425, L33 
Watson, D., Hjorth, J., Jakobsson, P., Xu, D., Fynbo, J. P. U., Sollerman, J., 

Thone, C. C, & Pedersen, K. 2006, A&A, 454, L123 
Yun, M. S., & Carilli, C. L. 2002, ApJ, 568, 88 
Zhang, B., & Meszaros, P. 2001, ApJ, 552, L35 



TABLE 1 
JVLA Observations 



Date 


St 


Frequency 


F„ 


(UT) 


(d) 


(GHz) 




2012 Aug 4.97 


0.93 


4.9 


44±8 






6.7 


25±7 


2012 Aug 5.99 


1.95 


4.9 


40±8 




6.7 


26±7 


2012 Aug 8.02 


3.98 


4.9 


45±9 






6.7 


25 ±8 


2012 Aug 14.94 


10.90 


4.9 


45 ±10 






6.7 


24±9 


2012 Sep 11.01 


37.97 


4.9 


<40" 




6.7 


<36" 



NOTE. — ° Limits are 3cr. 
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TABLE 2 

Optical and Near-IR Observations of the Host and Nearby Galaxy 



Filter 




Host Galaxy 


Nearby Galaxy 






AB Magnitude" 


f « 


AB Magnitude" 


Fv" 




(tan) 








d*Jy) 


;■ 


0.630 


25.9 ±0.2* 


0.16±0.04 i ' 


26.0 ±0.2 


0.15 ±0.03 


(' 


0.779 


24.8±0.15 


0.44 ±0.07 


24.4 ±0.1 


0.63 ±0.06 


Y 


1.021 


23.7 ±0.3 


1.20 ±0.40 


23.0 ±0.2 


2.30 ±0.46 


J 


1.235 


23.0 ±0.2 


2.19 ±0.44 


22.4±0.15 


3.98 ±0.60 


K s 


2.146 


22.0 ±0.1 


5.75 ±0.58 


21.8 ±0.1 


6.85 ±0.69 


Note. 


a ' 


These magnitudes and flux densities have been corrected for Galactic 



extinction of E(B-V) ~ 0.204 mag iSchlafly & Finkbeinei 2011). * This value has 
been corrected for an expected afterglow contribution of ~ 0.03 fijy. 
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• 



2012 August 4.27 





1 



2012 August 7.30 






- Subtraction 





































FIG . 1 . — Gemini-North ('-band images centered on the location of the Swift/XKT error circle (black circle; 1 .4" radius 90% containment) obtained at 5; S3 0.23 
d (left) and St fa 3.26 d (center), along with a subtraction of the two epochs (right). The subtraction reveals a fading optical afterglow (cross-hairs). The Chandra 
X-ray afterglow position coincides with the optical position (red circle; 0.54" radius 95% containment). Also seen is a galaxy located at a projected angular 
distance of S3 1.4" south-east of the host galaxy position. Each image is 20" on a side. 
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Time since burst (d) 

FIG. 2. — X-ray light curve, optical i'-band detection, and radio 5.8 GHz upper limits for the afterglow of GRB 120804A. The solid lines are an afterglow model 
fit using the formulation of Granot & Sari 1 2002), with E^ Im ~ 8 X 10 51 erg, n ss 10~ 3 cm" 3 , e e ss 0.3, eg Ri 0.1, and p fa 2.2. The low observed flux density in 
the optical band is suggestive of extinction, and can be explained with Ajj? st 2 mag (dashed red line). The dotted blue line indicates the predicted radio light 
curve for a model with fixed parameters of n = 1 cm" 3 and iiK.iso = 6 X 10 51 erg (i.e., matched to E^jeoi, that produces an indistinguishable fit in the X-rays (this 
model requires e e 0.25, eg 0.1, and p 2.1). Clearly, such a high density can be ruled out. The arrows mark the Chandra and XMM-Newton observations, 
including the time of a late Chandra observation. 
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Observed Wavelength fam) Observed Wavelength (^m) 



FIG. 3. — Left: Spectral energy distribution of the host galaxy of GRB 120804A (black circles: optical, near-IR, and radio) compared to a scaled SED of 
Arp220 (blue line) and an elliptical galaxy with a 2 Gyr stellar population (red line) at z = 1.3, Both SEDs are scaled to match the optical/near-IR photometry 
of the host galaxy. The inset shows a zoom-in on the optical/near-IR range, highlighting that both SEDs provide a reasonable fit in this wavelength range (with 
a somewhat better fit in ^,-band for the Arp220 SED). The gray symbols mark the photometry of the nearby galaxy, along with an elliptical galaxy model at 
z = 1.3. In both cases, the apparent steepening between the i- and F-band filters points to a similar redshift. Right: Comparison to the scaled SEDs of radio-quiet 
(red line) and radio-loud (red line) AGN at z= 1.3, as well as the SED of an elliptical galaxy with a 2 Gyr stellar population (green line). The AGN SEDs are 
scaled to match the observed radio flux density of the host galaxy. The radio-quiet AGN scenario over-predicts the optical, near-IR, and X-ray fluxes. On the 
other hand, the radio-loud scenario matches the radio flux density without violating the other measurements; the optical/near-IR emission is dominated by stellar 
emission. 



